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INTERACTIVE  SCENARIO  COMPUTER  MODEL  FOR  DOSFJ  RATES 
TO  AIRCREWS  IN  FLIGHT  THROUGH  NUCLEAR  DEBRIS  CLOUDS 


INTRODUCTION 

The  dynamics  of  radioactive  debris  fallout  from  a  nuclear  weapons  .burst 
has  been  extensively  modeled,  and  there  are  excellent  codes  to  predict  the 
distribution  on  the  ground.  These  codes  have  been  under  development  for  the 
past  20-25  years.  The  WSEG-10  (1)  and  DELFIC  (2)  are  notable  models  for  this 
purpose.  The  DELFIC  code  based  on  numerical  computation  is  the  most 
extensively  developed  dynamic  modex  and  is  used  as  the  benchmark  for  all 
analytical  work.  A  scenario  cf  comparable  strategic  importance  concerns  the 
calculation  of  dose  rates  to  aircraft  crewmembers  flying  through  nuclear 
debris  clouds  raised  aloft  by  ground  detonations.  This  scenario  would  occur, 
for  example,  in  Airborne  Command  Post  flights  of  the  Strategic  Air  Command 
immediately  following  a  nuclear  attack. 

The  penetration  of  a  nuclear  dust  cloud  is  hazardous  to  the  aircraft  and 
its  crewmembers.  Nuclear  debris  dust  is  also  abrasive  and  destructive  to  jet 
engines.  Gamma  radiation  from  the  dust  will  injure  personnel  as  well  as  the 
avionics  equipment.  This  report  concerns  an  interactive  computer  model  to 
assess  the  extent  of  the  radiation  hazard. 

Early  efforts  by  Patrick  et  al.  (3)  to  model  an  airborne  exposure 
scenario  were  based  on  empirical  approximations  to  the  nuclear  detonation 
cloud  dimensions  and  concentration  parameters.  This  empirical  model, 
however,  could  not  account  for  wind  parameters,  other  local  atmospheric 
conditions,  and  of  most  concern,  inhomogeneity  of  the  cloud  and  the 
fracticned  distribution  of  radioactivity.  These  parameters  have  a 
significant  effect  on  the  expected  dose  rates.  A  more  recent  model  by 
Hickman  (4)  is  based  on  WSEG-10,  a  revised  fallout  prediction  model  developed 
at  the  U.S.  Air  Force  Institute  of  Technology.  The  WSEG-10  model  developed  by 
Bridgeman  and  Bigelow  (5)  takes  into  account  fractionation,  wind  parameters, 
and  some  flight  parameters.  Yet  it  is  still  desirable  to  apply  the  full 
advantages  of  the  DELFIC  model  in  fly- through  scenarios.  Having  a  more 
detailed  structure  of  the  nuclear  dust  distribution  in  the  atmosphere  allows 
for  the  analysis  of  evasive  flight  trajectories,  and  loiter  scenarios  where 
the  aircraft  is  scheduled  to  fly  orbits  in  cloud  fields.  The  user  is  also 
assured  of  the  most  accurate  accounting  of  the  ground  zero  (GZ)  parameters, 
radioactive  dust,  winds  aloft,  and  nuclear  weapon  parameters. 

Implementation  of  the  full  capability  of  the  DELFIC  system  has  not 
heretofore  been  attempted  for  modeling  the  exposure  of  aircrewmembers  In 
fly-through  scenarios.  The  primary  criticism  is  that  the  program  i£  too 
expensive  to  run  and  that  the  model  is  difficult  to  learn.  In  spite  of  these 
circumstances  we  report  the  development  for  the  first  time  of  a  completely 
interactive  user-oriented  model  which  has  for  its  basis  the  benchmark  DELFIC 
code.  This  user-oriented  model  has  been  achieved  through  an  extensive 
adaptation  of  a  related  DELFIC  based  code  CASSANDRA  (6)  to  a  time-shared 
Digital  Equipment  Corporation  (DEC)  model  VAX  11/780  computer  system.  The 


adaptation  consisted  in  developing  a  large  entry  data  array  using  the  DEC 
Forms  Management  System  (FMS)  utility  and  in  developing  an  appropriate 
nuclear  scenario  modeling  theory.  A  complete  interactive  modeling  session  by 
an  individual  with  any  degree  of  computer  experience  can  be  accomplished  in 
less  than  20  min,  assuming  usual  computer  loads.  The  new  facility  can  be  a 
very  useful  aid  to  high  level  planners  and  modelers  of  nuclear  encounter 
scenarios.  Highly  detailed  analyses  can  be  rapidly  and  efficiently  performed. 
Studies  of  the  effects  of  variation  of  any  of  a  large  number  of  parameters 
are  possible.  An  example  of  the  list  of  parameters  that  can  be  easily 
specified  is  given  in  Table  1 . 


TABLE  1 .  PARTIAL  LIST  OF  PARAMETERS 


No.  Parameters 


1 . 

Height  of  burst  above  GZ 

2. 

Yield,  fission  fraction 

' 

3. 

Soil  type 

H. 

Particle  distribution  function 

5. 

Particle  solidification  temperature 

6. 

Wind  layers 

7. 

Atmosphere  variables,  humidity,  pressure, 

etc. 

8. 

GZ  altitude 

9. 

Any  number  of  aircraft  coordinates  at  any 
relative  to  GZ  and  to 

point  in  space 

10. 

Aircraft  air  mass  exchange  ratio 

11  . 

Cabin  dimensions 

The  following  sections  describe  the  theory,  program  structures  and 
procedures,  and  some  numerical  examples  of  the  interactive  model. 


THEORY 

The  nuclear  scenario  is  illustrated  in  Figure  I.  An  aircraft  i3 
designated  to  be  located  at  V^x*  ,yi,zi)  and  time  ( t i )  relative  to  GZ  and  to. 
Note  that  this  takes  into  account  the  aircraft  velocity,  and  that  time  to  is 
when  the  detonation  occurs.  The  critical  required  information  is  that  of  the 
dust  concentration  at  each  of  the  vertices  and  points  in  between  In  the 
flight  trajectory. 

The,  extraction  of  dust  concentration  computations  from  the  DELFIC  source 
code  is  available  in  the  CASSANDRA  code  (Cloud  Snapshots  of  Dust  Raised 
Aloft).  This  code  developed  by  Showers  and  Crisco  (6)  at  the  US  Army 
Armament  Research  and  Development  Command,  Aberdeen  Proving  Ground,  MD,  is 
designed  to  use  the  detailed  cloud  rise  module  of  the  DELFIC  system  which 
provides  cloud  dynamic  parameters  from  the  initial  time  to  stabilization.  The 
CASSANDRA  code  also  uses  the  DELFIC  atmosphere  model  and  particle  fall 
procedures.  The  stabilized  cloud  is  modeled  as  a  stack  of  wafers  of  given 
particle  sizes.  These  wafers  are  in  turn  transported  by  winds  and  deposited 
in  agreement  with  DELFIC  fall  dynamics.  The  CASSANDRA  code  was  developed  to 
assess  the  threat  caused  by  dust  to  reentry  vehicles  (RV)  reentering  the 


atmosphere  or  antiballistio  missies  (ABM)  departing  the  environment  of  a 
nuclear  detonation.  The  code  ideally  calculates  du3t  concentrations  for  a 
given  set  of  coordinates.  It  also  calculates  a  user  specified  2-D  map  of 
dust  in  a  plane  cutting  through  the  cloud.  With  the  availability  of  these 
calculations,  it  was  then  clear  that  fly-through  scenario  dose  calculations 
could  proceed. 

To  calculate  the  irradiation  dose  required  a  scaling  factor  related  to 
dust  mass,  the  assumption  of  two  radiation  sources,  the  outside  air 
(immersion  dose  rate),  and  internal  dose  due  to  dust  ingested  into  the 
aircraft.  Mathematically  it  is  assumed  that  the  aircraft  is  at  a  position 
vector  _tq  as  shown  in  Figure  2  and  that  the  irradiation  due  to  cloud  sources 
is  locally  homogeneous.  The  total  dose  rate  is  composed  of  two  factors, 
i.e., 


Total  dose  rate  ■  immersion  dose  rate  ♦  onboard  dose  rate 

The  immersion  dose  rate  (DOSER)  at  a  point  (r0)  and  time  (t)  is  given  by 

too 
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where  the  Way-Wlgner  (7)  time  dependence  i3  assumed,  t  is  in  hours  after 
detonation,  and  K1  is  a  factor  associated  with  the  energy  deposition  of  1  MeV 
gamma  ray  photons,  K,  -  4.88  x  10"10  rad  cm2/photon;  see  reference  8.  The 
summation  is  over  the  particle  size  class  J,  assuming  there  are  100  values. 
K(J)  is  the  radioactivity  associated  with  particle  class  J.  The  dust  density 
p(£0»  J)  a^30  depends  on  the  particle  class  J,  location  in  the  cloud  space 
£q  and  time  t.  For  the  purpose  of  computing  the  integral,  the  dust  density 
is  assumed  uniform  in  the  vicinity  of  the  integrated  region.  Thus, 
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Assuming  mathematically  that  the  integral  goes  over  Infinite  volume  results 
in  the  integral  having  the  value  in  spherical  coordinates  R  *  £  -  £0  * 


K3  is  the  atmospheric  absorption,  6.767  x  10“3  m"1  (*0.  Thus, 
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The  K(J)  parameter  Is  developed  no  follows.  Assume  tnat  the  total 
radioactivity  raised  aloft  is  530-gamma  MCi/Mt  of  fission  yield  and  Is 
normally  distributed  over  dust  particle  maos  (size)  (9).  This  takes  into 
account  fractionation  aa  discua3od  by  Bridgeman  and  Bigelow  (5).  Since  the 
radioactivity  1$  distributed  with  the  total  maos,  and  each  of  the  100 
particle  size  claaooo  carries  equal  mass,  then  the  dose  rate  becomes 

/oo 
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where  A  io  the  total  radioactivity,  A  -  7.06  x  1 0^2  y  .  g  (photori3/hour) 
where  Y  Is  the  yield  in  kilotons,  F  is  the  fission  fraction,  and  M  is  the. 
total-maoo  raised  aloft  (in  grams).  All  those  parameters  are  part  of  the 
DEL FIG  auction  of  the  model.  Given  a  sufficient  number  of  calculated  points, 
then  tho  total  dose  due  to  Immersion  can  be  computed  by  numerical  integration 
of  the  dose  rate  along  the  flight  trajectory  given  in  Figure  2, 

The  loci  of  tho  points  at  which  duct  concentrations  are  to  bo  computed 
arc  obtained  in  the  present  model  by  an  equal  space  incrementation  of  the 
Euclidean  distance  between  two  uaor-apecified  points,  l.e., 


n  ■  r.  a,  t  .i.  (n,  ~  n  )  ;  c  *  /,  Z,  3 


where 
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and  6  Is  an  incrementation  parameter;  from  experience  on  dust  concentration 
variations,  '•  -  '1,000  m  io  a  good  valuo.  The  time  nt  the  loci  is  similarity 
O'-vnpucud.  The;  total  dose  io  computed  by  numerically  integrating  the  average 
dose  rate  between  loci  and  the  transit  time  in  the  locuo-to-locuo  interval. 
It  remains  for  the  CASSANDRA  section  of  the  code  to  compute  the  100  particle 
dust  densities  at  each  intermediate  vertex  point. 

The  done  rate  due  to  the  dust  brought  onboard  io  computed  in  a  similar 
manner.  In  this  case  the  dust  brought  onboard  is  baaed  on  a  user  specified 
aircraft  parameter,  the  total  air  mass  flow  rate,  and  the  outside  dust 
donoity.  The  cabin  preoourlzirig  system  compresses  a  given  volume  of. outside 
ijir  and  with  it  carrying  a  certain  amount  of  radioactive  debris.  The  air 
volume  compressed  and  brought  in  (AVJNR )  is  given  by 

AV/A/R  e.  7Se>fF*)/  (9) 

where  FR  la  tho  maos  flow  rate  in  pounds  per  minute  and  ADZ  is  the  air 
density  at  Z  altitude  In  grams  per  cubio  centimeter.  The  ADZ  is  based  on  a 
polynomial  fit  to  atmospheric  data  (10),  The  total  dust  maos  (DBB)  is  equal 
to  the  product  of  tho  AVINR  parameter  with  the  duat  concentration  summed  over 


a 


particle  size  class  and  the  total  time  of  flight  between  the  Initial  locus 
and  the  point  in  question,  l.e., 
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With  the  total  dust  brought  onboard  thus  computed,  the  dose  rate  to  an 
aircrewmember  is  assumed  to  come  from  two  sources:  1)  the  dust  distributed 
throughout  tho  cabin  air  volume,  and  2)  the  dose  rate  from  a  potential  air 
filter.  The  cabin  is  modeled  as  a  right  circular  cylinder  of  (CL)  centimeter 
length  and  (CR)  centimeter  radius.  The  distance  of  the  aircrewmember  from  a 
filter  is  D  centimeter}  this  is  illustrated  in  Figure  3.  To  first 
approximation  the  model  assumes  a  fraction  (FR)  of  the  dust  goes  into  the 
cabin  air  and  the  remainder  stay3  in  the  filter.  The  cabin  air  dust  doso 
rate  (COUSTDR)  at  (i)  is  given  by 
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where  CV  »  ir  (CR)2cL  (cabin  volume).  In  cylindrical  coordinates  and  making 
the  approximation  exp[-kr]  =  1  over  the  dimensions  of  the  aircraft  gives 
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The-'f liter  du3t  dose  rate  (FDUSTDR)  is  simply  the  Inverse  distance  dependence 
of  the  remaining  dust  at  a  distance  D,  i.e., 

FAoSfA^lt)  *  ~  F*)  li. _ (Hd-  088  ft' )  m 
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These  dust  dose  •"'•tea  are  integrated  by  taking  the  mean  value  of  the  rate 
between  the  two  vertex  points  and  multiplying  by  the  time  interval.  In  the 
computer  program  the  air-  density  at  altitude  Z ( 1 )  is  computed  from  a 
polynomial  fit  to  a  standard  atmosphere  table.  The  computations  involved 
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hero  and  all  of  the.  Input  data  requirements  are  all  prompted  by  Interactive 
subroutines  aa  discussed  in  the  next  sootlon. 

Program  Structure  and  Procedure 

The  system  of  computor  subroutines  and  programs  which  comprise  thlo 
Interactive  computer  package  13  illustrated  in  Figure  The  original 

CASSANDRA  nomenclature  Is  maintained  for  Its  respective  programs}  all  other 
unita  arc  associated  uniquely  with  this  package  which  we  term  CASSIE. 
Extensive  reporting  features  arc  introduced  so  that  a  file  of  all  Input  data 
arrays  is  readily  available.  A  set  of  these  directories  (Tables  2,  3,  and  M 
comprises  the  total  system. 


TABLE  2.  DIRECTORY  [CASSIE.  SYSTEM] 


BOMB.COM; 1 

D0SER.EXE; 1 

CASSIE. OUT;  1 

CLOUD. OUT ;1 

POSER . DAT ;  1 

D03ERDAT.EXE; 1 

DUST.DAT ; 1 

M3DAT.EXE; 1 

DOSE R. OUT; 1 

D0SER0.DAT; 1 

M3PUST.FLB; 1 

MPDQ.DAT  j  1 

D03ERN.DAT  ;1 

M3DUST.EXE; 1 

MPDAT.EXE; 1 

MPDQ.OUT ;  1 

M3DUST.DAT; 1 

MEANDIAM.DAT; 1 

MPDQ.LAB; 1 

MPDQO.LAB;  1 

M3DUSTQ .DAT; 1 

MPDQ.FLB; 1 

MPDQ0.DAT; 1 

NEWVER.DAT;  1 

MPDQ.EXE; 1 

MPDQNEW.LAB;  1 

MXEC.OUT; 1 

MPDQNEW.DAT; 1 

MXEC.EXE; 1 

MXKC.DAT ; 1 

VERTEXO.DAT ; 1 

VERTEX. DAT; 1 

TABLE  3.  DIRECTORY  [CASSIE .SOURCES]  •; 

MPDAT.FOR; 1 

M3DUST .FOR;  1 

MPDQ.FOR; 1 

DOSERDAT .FOR ; 1 

DOSER.FCP.il 

MXEC.FOR ;  1 

M3PAT .FOR;  1 

TABLE  U.  DIRECTORY  [CASSIE. MISC] 

0  AS.BIHSUB.EXE;  1  CASSIE  .SUBJ  .OBJ ;  2 

OASSIESUB.FOR; 1  EDTINI .EDT ; 1 


The  execution  of  all  programs  in  this  system  (CASSIE)  is  controlled  by  a 
command  procedure  called  BOMB.COM.  It  is  started  by  typing  "@B0MB"  in 
response  to  the  operating  system  prompt.  Once  the  command  procedure  is 
entered,  tho  user  will  no  longer  need  to  worry  about  any  type  of  computer' 
commando,  and  all  keyboard  entries  will  be  in  the  form  of  answering  questions 
or  data  entry.  Upon  starting  tho  command  procedure,  the  soreen  will  clear 
and  the  menu  given  in  Figure  5  will  be  displayed.  The  user  will  then  be 
asked  to  enter  a  choice.  If  an  Invalid  choice  is  entered,  the  user  will  be 
told  and  reprotnpted  for  a  valid  response. 


MAIN  MENU 


You  have  the  following  options  from  which  to  choose: 


1) 

Create  the 

input  data 

for 

M3 DUCT 

2) 

Run  M3DU5T 

3) 

Create  the 

Input  data 

for 

MPDQ 

4) 

Run  MPDQ 

b) 

Create  the 

input  data 

for 

DOSER 

6) 

Compute  the  do3e  rates 

7) 

Exit  this 

procedure 

Figure  'j.  Cassie  nuclear  scenario  menu  display. 


When  the  chosen  program  terminates  normally,  the  message  "Normal  program 
termination."  will  appear  on  the  terminal.  After  a  normal  program 
'  'nation,  the  message  "Please  wait..."  will  appear  on  the  terminal.  During 

brief  pause  that  follows,  the  command  procedure  will  be  doing  some 

-al  housekeeping  chores.  These  include:  (1)  Renaming  old  versions  of 
files  that  were  updated  during  program  execution.  This  i3  accomplished  to 
restore  the  old  version  of  the  file  in  the  event  that  the  new  version  of  the 
fUe  is  incorrect.  (2)  Purging  any  excess  versions  of  files.  When  the 
housekeeping  is  completed,  the  menu  will  be  redisplayed  so  that  the  user  may 
enter  the  next  option.  Options  (1),  (3).  (5),  and  (6)  are  unions  to  the 
CASS IE  system  and  are  described: 

Options  1  and  3:  Creating  input  data  for  the  cloud  rise  and  dust 
concentration  modules 

Data  creation  for  these  modules  is  accomplished  through  the  use  of  the 
Digital  Equipment  Corporation's  FMS,  Through  the  use  of  this  utility,  a  user 
can  easily  enter  or  modify  program  data.  The  FMS  allows  the  user  to  enter 
data  via  a  computer  "form".  The  following  general  information  applies  to  all 
aspects  of  data  entry  during  terminal  session: 

1.  Data  fields  are  delimited  with  square  brackets  ("["  and 

2.  Pressing  the  (RETURN)  key  will  move  the  cursor  forward  to  the  next 
field. 

3.  Pressing  the  (BACKSPACE)  key  will  move  the  cursor  to  the  previous 
field. 

ii.  Pressing  the  (PF2)  key  on  the  VT100  terminal  will  provide  a  help 


Message  for  the  fiolu  that  the  cursor  is  on.  Pressing  the  (PF2)  key  again 
will  proviso  additional  help  information,  if  available. 

5.  Pressing  the  (PF4)  key  on  the  VT100  terminal  while  running  the  MPDAT 
program  will  produce  a  hard  copy  printout  of  the  ''form'*  a3  It  appears  on  the 
terminal  screen,  data  included. 

A.  Pressing  the  (LINEFEED)  key  will  delete  all  data  in  a  field. 

7.  Some  fields  are  validated  after  entry  to  insure  that  a  valid,  response 
is  entered.  In  the  event  that  an  entry  i3  invalid,  a  message  will  appear  and 
"the  user  will  need  to  reenter  a  value  for  that  field. 

The  data  entry  programs  were  designed  to  be  easy  to  use  for  a  user  of  any 
level  of  computer  expertise.  Help  messages  are  provided  to  help  the  user  to 
determine  what  is  being  requested.  Information  related  to  the  help  messages 
for  use  in  setting  up  Options  1  and  3  i3  derived  from  the  sample  data  cards 
of  Reference  6. 


Data  Flow  In  The  CASSIE  System 

The  interrelationship  of  files  and  programs  in  the  CASSIE  system  is 
given  in  Figure  4.  Data  for  the  M3DUST  program  is  created  through  the  use  of 
the  data  entry  program  M3DAT.  Input3  to  this  program  are  the  existing  M3DUST 
data  file  (M3DUST.DAT)  and  user  inputs  from  the  terminal.  The  library  that 
contains  the  computer  "forms"  used  for  data  entry  is  M3DUST.FLB.  Output 
from  the  data  creation  program  is  M3DUST.DAT,  an  updated  version  of  the 
original  data  file. 

Data  for  the  MPDQ  program  is  created  through  the  use  of  the  data  entry 
program  MPDAT.1  Inputs  to  this  program  are  the  existing  MPDQ  data  file 
(MPDQ.DAT) ,  a  work  file  (VERTEX.DAT),  and  user  inputs  from  the  terminal.  The 
library  that  contains  the  computer  "forms"  used  for  data  entry  is  MPDQ.FLB. 
The  VERTEX.DAT  work  file  contains  user-specified  vertices  that  were  used  to 
create  the  current  input  data  file  for  MPDQ.  Outputs  from  the  data  creation 
program  are  an  updated  VERTEX.DAT,  an  updated  MPDQ.DAT,  and  a  descriptive 
file  of  the  newly  created  Inputted  data.  This  file  is  named  MPDQ. LAB  and 
contains  the  input  data  along  with  a  description  of  the  data  item.  It  can  be 
printed  out  by  using  the  appropriate  print  command. 

Program  M3DUST  uses  as  input  the  data  file  M3DUST.DAT  and  creates  as 
output  MEANDIAM.DAT,  CASSIE. OUT,  and  CLOUD. OUT.  CASSIE. OUT  is  the  report 
file  from  M3DUST  and  can  be  printed  out.  CLOUD. OUT  is  a  binary  file  and 
CANNOT  be  printed.  Its  purpose  is  to  serve  as  input  to  program  MPDQ,  and  it 
contains  the  DELFIC  cloud  rise  information  for  the  computed  cloud. 

The  input  data  files  for  program  MDPQ  are  CLOUD. OUT  (created  by  M3DUST) 
and  the  MPDQ  data  file  MPDQ.DAT.  Outputs  are  DUST.DAT  and  MPDQ. OUT.  MPDQ. OUT 
■can  be  printed  on  the  system  printer  to  obtain  a  hard  copy  of  the  MPDQ 
report. 

Options  5  and  6:  Calculating  the  dose  rates 

Dose  rates  are  computed  through  the  use  of  program  DOSER.  This  program 
U3C3  as  input  DUST.DAT,  CLOUD. OUT,  VERTEX.DAT,  and  DOSERDAT.OUT.  VERTEX.DAT 
contains  the  user-specified  input  vertices.  DUST.DAT  is  created  when  the 
MPDQ  variable  COMBIN  is  equal  to  1.  This  file  contains  the  dust  amounts  for 
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a  given  point  in  3pace.  The  l'ile  la  separated  into  sections,  and  t no  number 
of  sections  is  dependent  upon  the  number  of  vertex  points  used.  Within  each 
section,  the  data  points  ropreaent  the  dust  amounts  for  each  particle  size 
class,  with  the  first  point  in  the  section  corresponding  to  the  first 
particle  size  class,  and  the  last  point  corresponding  to  the  last  particle 
size  class.  Option  5  calls  program  DOSERDAT  which  requires  user  inputs 
related  to  aircraft  characteristics.  OutDuts  from  DOSER  is  finally  the 
report  file  DOSER. OUT  which  can  be  printed  to  receive  the  hard-copy  report. 
This  is  the  final  output  product  of  the  interactive  session. 

CASSIESUBS  is  a  file  of  subroutines  used  by  both  MPDQ  and  M3DU5T.  The 
file  of  source  codes  for  the  subroutines  can  be  found  in  [CASSIE. SOURCES] 
while  the  compiled  subroutines  are  kept  in  [CASSIE. MISC]  (see  Tables  2,  3. 
and  4).  The  computer  "forms"  U3ed  by  M3DAT  and  MPDAT  to  input  data  are 
contained  in  the  form  libraries  M3DUST.FLB  (M3DAT  forms)  and  MPDQ.FLB  CMP DAT 
f  orms) . 


Numerical  Example 

The  model  can  be  easily  exercised  to  examine  a  scenario  associated  with  a 
1-Mt  yield  over  sandy  soil  conditions  mimicking  the  Dial  Pack  ground  zero 
conditions  and  atmosphere  (6).  For  the  example  treated  here  the  aircraft  is 
assumed  to  fly  through  at  an  altitude  of  13.500  m  with  two  stipulated 
vectors,  40  km  west  of  GZ  at  t+  -  784  s  and  40  km  east  of  GZ  at  t+  =,1234  s. 
Figure  6  shows  the  MPDQ  west-east  dust  concentration  looking  north  through 
the  X-Z  plane.  An  aircraft  figure  is  drawn  to  indicate  the  relative  position 
'a'"'  t+  =  1031  s.  The  dose  rates  were  computed  using  typical  parameters  for  a 
XC-135  aircraft.  A  printout  of  the  DOSER. OUT  file  which  contains  this  dose 
rate  data  is  giver,  in  Figure  7.  Note  that  the  immersion  dose  rate  as  the 
aircraft  leaves  the  cloud  drops  off,  whereas  the  onboard  radioactivity, 
although  small,  still  remains.  To  examine  the  significance  of  time  on  the 
dcse  rates,  Figure  8  gives  the  same  scenario  with  fly  through  at  t+  *  1  h  * 
3600  s.  Note  that  the  d03s  rates  have  fallen  off  sharply.  These  results 
appear  to  depart  somewhat  at  later  times  when  compared  to  previous  models 
(see  References  3  and  4). 

As  there  are  a  very  large  number  of  possible  scenario  parameters,  the 
ease  in  which  thi3  model  may  be  adapted  can  lead  to  numerous  analyses.  One 
example  not  discussed  is  the  use  of  GZ  wind  parameters.  Orbits  in  clouds  can 
bo  specified  by  prescribing  segmented  sections. 
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